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Evolution of primordial fluctuations in a Brans-Dicke type scalar-tensor gravity theory is com- 
prehensively investigated. The harmonic attractor model, in which the scalar field has its harmonic 
effective potential in the Einstein conformal frame and the theory relaxes toward Einstein gravity 
with time, is considered. The evolution of adiabatic initial perturbations in flat SCDM models is 
examined from the radiation-dominated epoch up to the present. We discuss how the scalar-tensor 
gravity affects the evolution of metric and matter perturbations, mainly focusing on the observa- 
tional consequences, i.e., the matter power spectrum and the power spectrum of cosmic microwave 
background temperature. We find that the early time deviation is characterized only by the large 
static gravitational constant while the late time behavior is qualitatively different from that in 
Einstein gravity because the time variation of the gravitational constant and its fluctuation have 
non-negligible effects. The attracting scalar-tensor gravity affects only small scale modes due to its 
attracting nature, the degree of which is far beyond the post-Newtonian deviation at the present 
epoch. 
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I. INTRODUCTION 

The existence of massless scalar partners associated 
with the tensor field of Einstein gravity is generically pre- 
dicted by the recent attempts toward unifying all elemen- 
tary forces in nature based on supergravity, superstrings 
0, and other higher dimensional gravity theories. In 
these theories, time variations of fundamental constants 
such as the gravitational constant and the fine structure 
constant are naturally introduced. A recent claim of the 
time varying fine structure constant from observations of 
QSO absorption lines M may be a piece of evidence. 

Scalar-tensor gravity theories, whose original version 
was proposed by Jordan || and Brans and Dicke 0] and 
were extended in a more general framework later Q , nat- 
urally provide coupling between the massless scalar fields 
and the tensor field of Einstein gravity. Scalar-tensor 
gravity theories are almost only possible alternatives to 
Einstein gravity. Moreover, scalar-tensor theories may 
supply a new approach to implementing the inflation- 
ary scenario (called extended inflation) ^-0. The time- 
varying gravitational "constant" which is their distinctive 
feature slows the inflationary expansion from exponential 
to power-law in time, and then the inflationary epoch 
has finite period, thereby solving the so-called "graceful 
exit" problem. Furthermore, scalar-tensor theories pro- 
vide a natural framework of realizing the time-variation 
of fundamental constants (gravitational constant) via the 



dynamics of the Brans-Dicke dilaton (for review see p0|). 

In the Jordan-Brans-Dicke theory [[§,0 (hereafter, we 
refer to it as the Brans-Dicke theory for simplicity) which 
is the simplest example of scalar-tensor theories, a con- 
stant parameter uj is introduced. In the limit u) — > oo, 
the gravitational constant can not change and Einstein 
gravity is recovered. Although scalar-tensor theories in- 
cluding the Brans-Dicke theory are compatible with Ein- 
stein gravity in several aspects, they have many devi- 
ations from it. Weak- field experimental tests in solar- 
system have constrained the post-Newtonian deviation 
from Einstein gravity, u> > 500 jllj , . Measurement of 
the signal time delay of millisecond pulsars or the light 
deflection of quasars may raise this limit ]l3| . 

In more general scalar-tensor theories, w can vary de- 
pending on the scalar field. In cosmological models based 
on such theories, it has been pointed out that there is 
generally an attractor mechanism that drives w to oo in 
the late cosmological epochs |14 . The nature of grav- 



* Electronic address: |nagata@th, nao.ac.jp 

tAlso at Max Planck Institute for Astrophysics, Karl- 
Schwarzschild-Str. 1 Postfach 1317 D-85741 Garching, Germany 



ity can be significantly different in the early universe. 
Hence, information on the different cosmological epochs 
may constrain such theories. A simple and natural ex- 
tension of the Brans-Dicke theory to the attractor model 
is the harmonic attractor model in which the scalar field 
has a quadratic effective potential of a positive curva- 
ture in the Einstein conformal frame. The analysis of 
big-bang nucleosynthesis (BBN) in this model jl5| re- 
stricts two parameters characterizing the potential (its 
curvature (3 and today's gradient oq). It is concluded 
that the BBN limit on the possible deviation from Ein- 
stein gravity (2u>q + 3 = cyq~ 2 ) is much stronger than 
the present observational limits in large /3(> 0.3) mod- 
els. Aside from BBN, we have another source of the 
information about the early universe that is the cosmic 
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microwave background (CMB). The trace of primordial 
fluctuation can be seen clearly in the CMB anisotropy 
spectrum where the information on the early universe up 
to the last scattering time (z ~ 1000) is integrated on 
the acoustic peaks. The precise CMB data which will be 
provided in near future , [ [l7j should further constrain 
the allowed parameter region. Toward probing gravity 
theories by CMB data, we examine the cosmological per- 
turbation evolution in the context of scalar-tensor gravity 
theories. 

Our primary goal is to clarify the growth process of 
perturbations in the scalar-tensor cosmology and under- 
stand the influence on the resultant observable power 
spectra. We also revisit the perturbations of the Brans- 
Dicke theory which has been investigated by several au- 
thors [|l8| [n| |o). Although the contemporary plausi- 
ble model is a vacuum dominated model, we adopt flat 
standard CDM (SCDM) models for three reasons. The 
first reason is that the universe necessarily experiences 
matter domination once at least. Therefore, we should 
begin by revealing the effect of the scalar-tensor grav- 
ity up to the matter-dominated epoch without confusing 
curvature or vacuum energy effect with it. The second 
reason is that the dependence of the scalar evolution on 
the matter density (floh 2 ) is quite simple as shown in 
Section IV. The final reason is that, if we take account of 
the current constraint for u>q , the attracting scalar-tensor 
gravity virtually becomes Einstein gravity before vacuum 
domination and hence any extra process originating from 
the scalar-tensor gravity would not occur after vacuum 
domination. 

This paper is organized as follows: In Section II, the 
field equations and the background cosmological evolu- 
tion equations in the scalar-tensor theory are described. 
In Section III, we present the cosmological perturbation 
theory based on the scalar-tensor gravity. We demon- 
strate the numerical calculations in Section IV. The re- 
sultant matter power spectra and the CMB temperature 
anisotropy spectra are illustrated in Section V. Finally, 
some conclusions are in Section VI. In Appendix, the an- 
alytic approximate solutions of the background and the 
perturbation equations in the Brans-Dicke theory are dis- 
cussed. 



A. Field equation 

The action describing a general massless tensor- 
monoscalar theory is 



S 



1 

16^ 



+S m [i/j,g^}, (1) 



where R — g^R^ v denotes the scalar curvature of the 
metric <? M „. The last term in Eq.(Q) denotes the action of 
the matter which is a functional of the matter variable tp 
and the metric g^ u . Lu((f>) is a function of 4> and it works 
as the coupling parameter of the nonminimal coupling. 
We redefine the scalar field so that it is dimensionless 



G' 



(2) 



where G is the Newtonian gravitational constant mea- 
sured today. The field equations are 
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{V M 0V^-- 5 ^(V0) 2 } (3) 



and 



R 



2-Dc 



V T v 



1 dhJ 



(V0) 2 , 



0. 



(4) 
(5) 



where V M denotes the covariant derivative defined by the 
g^, □ = g M,y V M V„, and T" is the matter stress energy 
tensor defined as follows 



2 ss m [ ^,g^\ 



(6) 
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The energy- momentum conservation is satisfied as Eq.(||) 
because ordinary matter fields do not directly couple to 
the scalar field. Combining Eqs.(j^)and(Q), the equation 
of motion for d> becomes 



1 {8*GT;-£(V*) a }. (7) 
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II. MODEL 



In this section, we describe the scalar-tensor theory 
and the (unperturbed) cosmological model based on it. 
We do not assume a specific form of u>(<fi) as a function of 
the scalar field (f> until performing numerical calculations. 



B. Background cosmological equation 

The unperturbed cosmological spacetime metric is 

ds 2 = a 2 (77)(-(V + HjdJda?), (8) 

where jij is the metric on the comoving homogeneous 
isotropic 3-space of a constant spatial curvature K. The 
components of the unperturbed stress energy tensor are 

To° = -P: I? = 0, T]=p 7 j, (9) 



3 



where p and p are the total energy density and pressure, 
respectively. The background evolution equations are 



P 
K 



-3-(p+p), 
a 
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a <j) 



Lll 



2u 



U 

l—^Ga^p-Zp) 



U3 (<t>^ 2 

6 



> duj 



(10) 
(11) 
},(12) 



where the prime denotes a derivative with respect to 
77. Although Eq.(|l^) is the same as that in Einstein 
gravity as noted previously, the expansion rate a' /a is 
given by Eq.(O). Non-relativistic matter always con- 
tributes to the first term in the right hand side of Eq. ( 12 ) . 
The requirement that today's gravitational constant be 
in agreement with Newton's constant determines the 
present value of <fi as 



4 + 2uj 
3 + 2uj ' 



(13) 



where <fio and luq denote the present value of <f> and ui (</>), 
respectively. 



III. SCALAR-TENSOR COSMOLOGICAL 
PERTURBATION 

Following the gauge invariant formulation of the 
cosmological perturbation theory in Einstein gravity 
pl[.p2|], we formulate the gauge invariant perturbation 
theory in the scalar-tensor gravity. 



A. Perturbed quantities 

Perturbations of scalar quantities can be expanded by 
the harmonic functions Y defined as 



(A + k 2 )Y = 0, 



(14) 



where —k 2 represent their eigen values of Laplace- 
Beltrami operator A on the comoving homogeneous 3- 
space. Scalar type components in perturbations of vector 
quantities are expanded by 

Yt = -k~%, (15) 

where u denotes a covariant derivative associated with 
7y . Those in perturbations of tensor quantities are ex- 
panded by 



and 



Y 

1 1.1 



(16) 



(17) 



which correspond to the trace and the traceless part, re- 
spectively. In this paper, we will omit the summation 
symbols and the eigenvalue indices because there is no 
coupling among the different fc-modes. Denoting the per- 
turbed metric by g^, the components of the perturbed 
metric are written as 



500 = -a 2 (l + 2AY), 
90] = -a 2 BYj, 

gn = a 2 {(l + 2H L Y) lij +2H T Y ij }. 
The perturbed scalar field is 

4>= + X Y. 



(18) 
(19) 
(20) 

(21) 



The perturbed energy momentum tensor Tj^ has the fol- 
lowing components: 



1 Q 

J 

fj 

rpi 



-p(l + SY), (22) 

{p + p){v-B)Y h (23) 

-{p + p)vY\ (24) 

p{(i + ^ L r)7] + 7r T r;}, (25) 



where S, v, ttl and ttt are the perturbations of the en- 
ergy density, the spatial velocity and the isotropic and 
the anisotropic pressure, respectively. All of the above 
perturbation quantities are functions only of time. When 
we consider the background and the perturbed universe, 
there is gauge freedom to fix the correspondence between 
the spacetime points in the two spacetimes. Denoting a 
scalar type infinitesimal gauge transformation by the vec- 
tor 



TY, LY j ), 



(26) 



where T and L are arbitrary functions of time which has 
the same order as perturbation variables, we obtain the 
changes of perturbations under the gauge transformation 
as 



.4 
B 
Hl 



A-T'- (a'/a)T, 
B + L' + kT, 
H L - (k/S)L- {a'/a)T, 
H T + kL, 
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v 

7TT 



and 



v + L', 

- (p'/p)T, 

TTT, 



= x-4>'t, 



(27) 
(28) 
(29) 
(30) 

(31) 
(32) 
(33) 
(34) 



(35) 



where the variables with bar represent the gauge trans- 
formed values. To treat the evolution of perturba- 
tions without unphysical gauge modes, we deal only 
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with gauge-invariant variables. We can construct gauge- 
invariant variables by combining perturbation variables. 
For one perturbation variable, there are several gauge 
invariants each of which corresponds to the perturbation 
in a specific gauge. Here, we introduce two gauge in- 
variants for each of the metric, the scalar field and the 
density perturbations. The one denoted by subscript s 
corresponds to the perturbation in the Newtonian shear 
free gauge (where Ht = B = 0). The other denoted by 
subscript c corresponds to the perturbation in the total 
matter comoving gauge (where B — v). Let us proceed to 
the definition. We can employ two independent variables 
^ and $ for metric perturbations. The gauge invariant 
perturbation variables associated with the (0-0) compo- 
nent of the metric are defined as 

* s = A + k- r {a'/a)(B - fc^iV) 

+k- 1 (B' -k^Hr"), (36) 

* c = A + k^ia' /a){B - v) + k^ 1 {B' - v'). (37) 

These correspond to the perturbations of the gravita- 
tional potential in the two gauges. On the other hand, 
the gauge invariants associated with the spatial compo- 
nent of the metric are 

* s = H L + ^H T + k-^a' /a){B ~ k^Hr 1 ), (38) 

$ c = HL + ^HT + k^ia'/a^B-v), (39) 

which represent the perturbations of the intrinsic spatial 
curvature in the two gauges. Similarly we introduce two 
gauge invariants of the scalar field perturbation 

X s = x ~ k~ l 4>' (k~ l Ht' — B), (40) 
X c = x-A~V(«-B). (41) 



The gauge invariants of the energy momentum tensor 
perturbations are defined as follows: 

A s = 5 + 3(1 + w)(a' /a)k- l {k~ l H T ' - B), (42) 
A c = 5 + 3(l + w)(a'/a)k- 1 (v - B), (43) 



V = v-k~ x H T ', (44) 



r = tt l - ^5, (45) 
w 



n = 7T T , (46) 

where 

w = p/p, (47) 

c\ = p'/p'. (48) 

The conventional photon temperature fluctuation vari- 
able is 

©0 = ^A S7 , (49) 

where the subscript 7 means the photon component. 
B. Perturbation equation 

The perturbation equations of Eq. (0) are 



3(-) * s _ 3-$; - (fc 2 - ZK )*, 
\ a J a 



8ttG 



d(j> a 2 4> 



-*.s 



(50) 
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ktu 
a 2 4> 



(51) 
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2^"* s + + 2-* s - 2$' s ) - - -X' s - ^-X, 

(a J a 6 



2a? 



' duj 



a 2 4> 



(52) 



k 2 
a 



f -i 8vrG fc 2 



(53) 



The perturbation equation of Eq. (0) is 



X'' + 2-X' + k 2 X s - 20"* s - <£'(*'.. + 4-4- s - 3$' ) 
a a 



2w + 3 



8ttGp{(1 - 3c 2 )A s - 3wT} - —X^ 2 ^ 2 - 2— a - 2 {4>' X' s - 



r* s )-2— a- 2 X s ((i>" + 2-0') .(54) 
dtp a 



r 



The perturbation equations of Eq. (||) are 

A' c - 3^jwA c = -3(1 + w)$c - (l+w)kV 



r' (-)v = w 

k 



fc r 9 A 2 / 3if , 



A- 2 



}. (56) 



From Eq. fl56| ) , the relation between the potential pertur- 
bation and pressure perturbations becomes 



Employing the total matter gauge variables and using 



These are the same as those in Einstein gravity. Of 
course, the Boltzmann equations for photons and neu- 
trinos are also not modified. Hence we omit them here. we obtain the alternative expression of Eq.(E3) that is 



J 



K + \K + k 2 X c - 20"* c - 4>'(% + £* c ) + kj>'(l ~^§)V+ ~{(^)Xc + 0'* c - x>} K \ .v, 



2lo + 3 



^ ■'•■^> " ' ' ~~^> 2 -2—a- 2 (0'X' c -<f> ,2 ^ c )-2—a- 2 X c (0" + 2-0') .(59) 

' dtp 



8ttG P {{1 - 3c 2 ) A c - 3wT} - -r-~X c a 



3w 

. a' 
a 



Combining Eqs.(|50|), ( |5l|) and (|53"|), we obtain the gen- follows: 
eralized Poisson equation in the scalar-tensor theory as 
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The gravitational redshift factor in the scalar-tensor the- 
ory is written as: 

$ s - * 5 = I> + T , (61) 



T^, and are the contribution of matter perturbations 
and the scalar field perturbation, respectively. If the 
background spacetime and the initial conditions for per- 
turbations are given, we can solve the above equations. 
The initial conditions for perturbations (solutions on a 
superhorizon scale in the radiation-dominated epoch) are 
discussed in the next section. 



IV. NUMERICAL CALCULATION 

Let us now numerically solve the evolution of the back- 
ground universe and perturbations. We apply our formu- 
lation to flat SCDM models. The matter components are 
baryons, cold dark matters, photons, and three species 
of massless neutrinos. (In cold dark matter dominated 
models, massive neutrinos do not significantly change the 
evolution of 4>-) in our numerical calculations, the func- 
tional form of us(4>) is 

2co + 3 = {a 2 - pln^/fa)}' 1 . (64) 

This u>(4>) corresponds to the harmonic effective poten- 
tial of the scalar field in the Einstein conformal frame. 
ao and (3 in Eq. (|64|) are today's potential gradient and 
curvature, respectively [|l4|] , jl5|] . If (3 = 0, this model is 
reduced to the Brans-Dicke theory. Moreover, the model 
with ao — > 0, [3 — is Einstein gravity. To realize the cos- 
mological attractor mechanism, we treat the parameter 
region (3 > in this paper. Although some parameter re- 
gion is already excluded by solar-system experiments and 
the BBN analysis, we often adopt unrealistically large 
values of ao and (3 to see the qualitative parameter de- 
pendence of the results. 



I 

where 



(62) 
(63) 

I 

A. Boundary condition 

First, we describe the boundary conditions for back- 
ground variables. We employ the conventional definitions 
of cosmological parameters 

n Q = Po /p c , n K = -K/a 2 Hl p c = 3H%/8ttG. (65) 

Unless otherwise noted, we fix today's baryon density 
(f2&(= Pbo/Pc)) and Hubble parameter (h) to 0.03 and 
0.7, respectively. (In this paper, we do not discuss the 
light element abundance synthesized via BBN in the 
scalar-tensor cosmology. The helium mass fraction is 
fixed at 23% in our calculations.) Then, the density of 
CDM is determined by f2 c o = £lo — fifco- In the scalar- 
tensor theory, the relation between the cosmological pa- 
rameters is as follows 

^ + n K -S+^Si = l, (66) 
0o 6 

where S = 0o/0o-ffo an d = \-§^<i>- We set Hk in the 
above equation to be zero since we treat only flat mod- 
els. As mentioned in Section II, </>o must satisfy Eq.([l3"|). 
In the scalar-tensor theory, Qq is not necessarily unity 
even in flat SCDM models. There exists the degree of 
freedom to determine the amplitude of the vacuum solu- 
tion mode of cf>' (oc a~ 2 (2w + 3)"^) at the present time. 
As in the precedent study |2(i[j , we assume that such a 
mode is negligible ever since the initial time of our cal- 
culations. Then, from the particular solution during the 
matter dominated epoch p4| , the relation that f2o and S 
must satisfy in < [3 < 3/8 models becomes 

_ (l+qg)(l-ayg/3) 2a d<p 
iZ °~ (1 + ao^o) 2 ' (1 + ao^o)' 1 ' 



2w + 3/0' 



2<P V <j) 



{x' + z^Xc-^} 
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where d(p = -3a (l - r)/{A0) and r = (1 - 8/3/3) 1 / 2 . 
In this paper, we set the initial time z ~ 10 8 so that the 
modes observed today as large scale structures are still 
before their horizon entry. Moreover, the energy of elec- 
trons and positrons already have no effect on <j> evolution 
at that time. We search for the numerical solution of <j> 
satisfying the above boundary conditions by iterating the 
initial value of <f), 4>ini ■ The initial value of <p' to a given 
4>i„i which is composed of the particular solution is 
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a eq \ 
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(68) 



Here we assume that <p (and ui (</>)) are almost constant 
and the <j)' terms in Eq.(ll) are negligible at the initial 
time. In fact, the factor a/a eq extremely damps its am- 
plitude and this is consistent with these assumptions. 

Next, the initial conditions for perturbation variables 
should be noted. Here, the initial condition means the so- 
lution on a superhorizon scale in the radiation-dominated 
epoch. To obtain the initial solution of X c , we employ 
the similar procedure as applied to </>'. In addition, for 
matter perturbations, we adopt the same adiabatic initial 
conditions as those in Einstein gravity. Then, the homo- 
geneous solution of Eq.(^9|) consists of a constant and a 
decaying mode, whose contribution to the gravitational 
potential perturbation decreases with time. Neglecting 
them and using the initial solutions for matter pertur- 
bations, we obtain the growing particular solution of X c 
that is 



X c = 



5 + 3/„ 



32c<ji. 



7 + 5/„ 



(69) 



where f v — p v jp 1 . This solution contains neutrino 
anisotropic stress contribution. Consistently, this initial 
solution of X c does not modify the initial solutions for 
matter perturbations because the contribution of X c (and 
(/)') to metric perturbations is at most comparable to that 
of CDM perturbation 0(a/a eq ). 



B. Background 

Here we compare the time evolution of background 
variables in some typical models. 

Fig.|l| shows the ag dependence of <j> evolution and Fig.^ 
shows the /3 dependence of it. (f> is frozen during the 
radiation-dominated epoch and begins to roll down at the 
matter-radiation equality time to realize the Newtonian 
gravitational constant at the present. As increasing a$ 
or /3, we obtain smaller initial <fi- I n the Brans-Dicke 
(constant lu) models, the late time slope is monotonic. 
On the other hand, in the attractor (variable to) models, 
the time variation of uj{4>) which regulates the amplitude 
of 4>' (proportional to (2uj + 3) _1 ) causes a steep slope 
immediately after the equality time and a gentle slope 
near the present. For reference, the flgh 2 dependence of <j> 
evolution is illustrated in Fig.[| The shift of equality time 



simply decreases the initial value of <f> without altering the 
late time evolution. Now, even if we employ a vacuum 
dominated model, the period from vacuum domination is 
quite short in units of In a and hence the only influence 
on cj> is to decrease in today's matter density. Thus, in 
a vacuum dominated model, the history of <j) is modified 
only before the equality time as the low density SCDM 
models in Fig|^. 

Fig.|| compares the time evolution of uj{(j)) in the at- 
tractor models. u> can vary only after the equality time 
because, even in variable w models, the variation of <p 
is suppressed before the equality time. The late time 
power-law behavior is well described by 



2cj + 3~ (2w + 3)f— ) 
\a J 



|(l-r) 



(70) 



r is defined in the previous section. 

Fig.|5| illustrates the time evolution of cf>' /{a 1 /a) in the 
attractor models which also represents the variation of <f> 
during a logarithmic interval of 1 + z. (It should be noted 
about ao dependence that the whole amplitude of these 
curves is nearly proportional to ao 2 -) The variability of 
d> is summarized as follows: 



'(I)"' 



d\na 



1 



RD 

MD 



(71) 



The initial slope represents the ratio of matter energy 
density to radiation energy density. Then, the motion 
of 4> is negligible and does not affect the spacetime met- 
ric through either of the minimal and the non-minimal 
coupling. Hence, the scalar-tensor gravity effectively be- 
haves like Einstein gravity and the only difference is the 
magnitude of the gravitational coupling " constant" . The 
larger gravitational constant results in the larger expan- 
sion rate and therefore the smaller horizon length. On 
the other hand, in the matter-dominated epoch, 4> can 
move in the absence of the suppression factor 0(a/a eq ). 
As a result, the time variation of 4> changes the nature 
of gravity qualitatively. For example, a is not propor- 
tional to rf . Especially in the attractor models, such 
effects appear only in the intermediate epochs between 
the equality and the present because the deviation from 
Einstein gravity is in proportion to uj{4>)~ 1 . Thus, the 
effect of scalar field dynamics is most significant around 
the equality time. 



C. Perturbation 

In this section, the evolution of perturbation variables 
is shown. When comparing the perturbation spectra, we 
adopt scale invariant models without discussing the gen- 
eration mechanism of perturbations. The amplitude of 
the curvature perturbation ($ s ) in each model is initially 
indistinguishable. Hence, a comparison between models 
shows directly the difference in their evolution. 
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FIG. 1: The time evolution of <j> in the Brans-Dicke models 
with SCDM parameters where Oo is determined by Eg. (|67|) . 
The discrepancy at the present is due to the different ui. 
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FIG. 3: The density dependence of <f> history in the attractor 
models with SCDM parameters, fio is determined by Eq.(^) 



and Q. b h 2 is fixed to 0.0147 
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FIG. 2: The time evolution of (j> in the attractor models with 
SCDM parameters where f2o is determined by Eg. ([i^) . Note 
that the variation of u makes an inflection point around the 
equality time (cf. Fig.^). The alignment at the present time 
is due to the equality of u>{<j>o). 



1. Scalar field perturbation 

Fig.|| and Fig.0 show the early and the late time be- 
havior of the scalar field perturbation, respectively. The 
analytic estimation is presented in Appendix. Before 
the horizon entry, the scalar field fluctuation cumulates 
around the dense region and X c always grows. Then, 
the amplitude is approximately proportional to the fac- 
tor which regulates the nonminimal interaction. 
While, after the horizon entry, the field fluctuation flows 
out and its growth is suppressed. Then, X c oscillates at 
the same frequency as radiation fluid in the radiation- 
dominated epoch or converges to some amplitude in the 
matter-dominated epoch. 

As in the background, the contribution of the scalar 



FIG. 4: The time evolution of ui(cj>) in the attractor models 
with SCDM parameters. 
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FIG. 5: The time evolution of 4>/H in the attractor models 
with SCDM parameters, where <j> — \-^d>- The maximums 
correspond to the inflection points in FigM 
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FIG. 6: The early evolution of X c on a scale k = 10 Mpc~ 
in the scalar-tensor models with SCDM parameters. Normal- 
ization is arbitrary. Their amplitude is proportional to lu^ 1 
(or LUini -1 ). The discrepancy of their horizon entry times is 
due to the difference in 4>i„uiai- 
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FIG. 7: The late time evolution of X c on a scale k — 0.003 
Mpc -1 in the scalar-tensor models with SCDM parameters. 
Normalization is arbitrary. These modes enter the horizon 
around z ~ 10 2 . Their amplitude is proportional to w(</>) -1 
and the convergence at the present represents the same u)q . 



field to the metric is initially negligible. Until the 
mode enters the horizon, it becomes larger and reaches 
0(ct>i n j -1 ) fraction around the equality time. It is shown 
in Fig.^J and Fig.|| where the feature of the largest scale 
mode is quite similar to the respective curve in FigJE). Im- 
mediately after the horizon entry, the contribution sud- 
denly drops and becomes negligible. This potential drop 
at the horizon entry is the characteristic event caused 
by the scalar field fluctuation especially in the matter- 
dominated epoch during which the horizon entry does 
not involve the potential decay in the Einstein model. 
Because of the late time convergence to Einstein grav- 
ity, large scale modes, as well as small scale modes which 
enter the horizon long before the equality time, are not 
affected by the scalar fluctuation in the attractor model. 
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2. Metric perturbation 

The metric evolution directly affects the evolution of 
matter perturbations. Its evolution can be seen through 
the observations of the present matter structures. For ex- 
ample, the turnover scale in the matter power spectrum 
corresponds to the horizon scale at the equality time. 
CMB observations provide more detailed information 
about the gravity evolution. Therefore, we pick up some 
scale modes which appear on the CMB anisotropy spec- 
trum as the characteristic structure, namely the modes 
corresponding to the acoustic peaks. The evolution of 
the first and the second acoustic peak modes at the de- 
coupling time is displayed in Fig.|l^ and Fig.[ll] where the 
large scale limit modes are also compared for reference. 
The initial behavior is model independent because the 
scalar-tensor gravity behaves like Einstein gravity. In the 
Brans-Dicke model, the late time evolution, except for 



FIG. 8: The evolution of -T^,/($ s - * s ) on some typi- 
cal scales in the Brans-Dicke model with SCDM parameters. 
Small scale modes are not affected by the scalar fluctuation. 



the horizon entry, is flat as discussed in Appendix. Also 
in the attractor model, the late time evolution gradually 
becomes flat because of the convergence to Einstein grav- 
ity. The large scale limit mode especially is not affected 
by the difference in the horizon entry and its late time 
deviation is determined by u>{4>)~ ■ After the equality, 
the decay widths of peak modes at a time generally be- 
come smaller in the scalar-tensor models than those in 
the Einstein model (see e.g. the second peak mode in 
Fig.|l^ or Fig.pl]). We point out two reasons for it. One 
reason is that the non-relativistic matter amount (Oo) 
is larger in the scalar-tensor models. The relative con- 
tribution of radiation to metric perturbations is smaller 
at the same z and hence the decay caused by radiation 
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FIG. 9: The evolution of — T^/(<I> S — ty s ) on some typical 
scales in the attractor model with SCDM parameters. There 
exists the scale which is most affected by the scalar fluctua- 
tion. 

pressure becomes smaller. The other reason is the mo- 
tion of <fi. As explained in Fig.|l2|, if <j> is constant, a peak 
mode in each model enters the horizon at the same time. 
However, since (f> moves with time after the equality, the 
entry time of the peak mode in the scalar-tensor model 
is shifted to lower z. This time lag also makes the decay 
width smaller. Furthermore, during the early phase of 
potential decay (approximately from the horizon entry 
up to the first maximum compression of photon fluid), 
the contribution of the scalar induced decay occurring 
coherently to the decay deriving from radiation pressure 
is not negligible. For instance, the potential of the first 
peak mode in Fig.|l^ or Fig.jll] catches up with that in 
the Einstein model around z ~ 10 3 . 

The late time evolution of the power spectra of the 
gravitational potential in the large scale in the scalar- 
tensor models is displayed in Fig.[l3] and Figjl^. Since 
that scale enters the horizon much later after the equality 
time, the scalar induced potential decay can be clearly 
seen. As noted above, the decay width is decreasing with 
time in the attractor model and hence the large scale 
modes are not affected significantly. 



3. Matter perturbation 

In the radiation-dominated epoch, the difference in 
the magnitude of the gravitational constant alone affects 
matter evolution. It is clearly seen in Fig.|l5| which ex- 
hibits the acoustic oscillation of photon fluid at z = 10 6 . 
As explained in Fig.|l2|, the peak mode in each model en- 
ters the horizon at the same time due to the constancy 
of the gravitational constant (</>). Therefore, the degree 
of the peak boost is equal for each model. Parallel trans- 
ports make these curves identical except in their diffusion 
cut off. After the gravitational constant begins to vary, 
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FIG. 10: The evolution of "^s — on some scales in the 
Einstein and the Brans-Dicke model with SCDM parameters. 
Normalization is arbitrary. 
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FIG. 11: The evolution of "I> s — ^ s on some scales in the 
Einstein and the attractor model with SCDM parameters. 
Normalization is arbitrary. 



the peak heights generally become lower than those in the 
Einstein model. Exceptionally, during the first acoustic 
compression, the mode experiences also the enhancement 
caused by the scalar induced potential decay. 

The observable oldest universe is the photon decou- 
pling epoch. Fig.[l 7 j and Fig.[l7| show the photon and the 
potential power spectra at z = 1000. It should be noted 
that z of the recombination time is not shifted because 
it is determined almost only by the photon tempera- 
ture. On very large scales (commonly called Sachs- Wolfe 
plateau), the fluctuation amplitude in the scalar-tensor 
models is slightly larger than that in the Einstein model. 
The deviation represents the value of cj" 1 at that time. 
The small scale potential height is, of course, higher in 
the scalar-tensor models, which implies the smaller hori- 
zon length at the equality time. On acoustic peak scales, 
in addition to peak locations which represent the hori- 
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FIG. 12: The schematic diagram illustrating how the varia- 
tion of <f) reduces the potential decay width, where ao/a = 
1 + z. For instance, let Xe and Ast be the second peak scales 
at some redshift in the Einstein and the scalar-tensor model, 
respectively. Due to the smaller <f> in the past, the horizon 
expansion in the scalar-tensor model behaves as lower solid 
line in the figure rather than dot-dashed line and hence the 
horizon entry is shifted to lower z than that in the Einstein 
model. 
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FIG. 13: The power spectra of gravitational potential $ s — 
tys at some redshifts in the Brans-Dicke model with SCDM 
parameters. Normalization is arbitrary. The potential decays 
to some extent as the mode enters the horizon. 

zon length at the time, we can see three features: the 
first Doppler peak height, the second peak height and 
the diffusion envelope. Let us begin with the diffusion 
scale. The diffusion length is also dependent on the hori- 
zon length. In fact, it is defined as: 

l»=(—f , (72) 

where Z#,n e , and gt are horizon length, free electron 
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FIG. 14: The power spectra of gravitational potential & s — ty 3 
at some redshifts in the attractor model with SCDM param- 
eters. Normalization is arbitrary. The decay width becomes 
smaller as to becomes larger. 



density, and the Thomson scattering cross-section, re- 
spectively. And a rec is the scale factor at the beginning 
of hydrogen recombination. The mean free path is model 
independent and hence the diffusion efficiency depends 
only on In 2 - Therefore, the higher diffusion envelope 
in the scalar-tensor models also derives from the smaller 
horizon length. On the other hand, single peak inside the 
envelope in the scalar-tensor models is lower than that in 

the Einstein model because the ratio of the peak scales 

_ i 

to the diffusion length is, in turn, proportional to Ih 2 ■ 
Aside from the diffusion effect, acoustic peak heights in 
the scalar-tensor models generally become lower after the 
equality time due to the smaller potential decay of the 
peak modes. This results in the lower second peak in 
Fig.jlj] and FigJlT]. This difference in the second peak 
height indeed originates from the difference in gravity, 
which can be confirmed by trying different Qb models. 
Concerning the first peak scale, as noted in the previ- 
ous section, the scalar induced potential decay has non 
negligible effect. It occurs coherently to ordinary poten- 
tial decay and makes the total decay width larger. This 
boosts the first Doppler peak higher because the mode is 
in its first maximum compression phase at the decoupling 
time. (Precisely, ordinary potential decay is triggered by 
a "sound" horizon entry. Thus the degree of this enhance- 
ment correlates with the sound velocity of photon-baryon 
fluid. In larger f^ft 2 models, the enhancement of the first 
peak becomes more significant.) 

Also after the decoupling time, the scalar induced 
potential decay affects matter evolution at the horizon 
entry. It is found that it makes a bump in the po- 
tential transfer function (Figjl^ and Fig.|l4|) and the 
photon root mean square fluctuation (|9o + ^s\rms = 

1 O + * s | 2 + 57+1 ) spectrum (Fig.fTJ) around the hori- 
zon scale of the time. In Fig.p~8L the amplitude on super 
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FIG. 15: The power spectra of effective temperature fluctua- 
tion at z — 10 6 in the Einstein and the scalar-tensor models 
with SCDM parameters. Normalization is arbitrary. Before 
4> begins to move, the only difference is that the gravitational 
"constant" in the scalar-tensor models is larger than Newton's 
constant. 
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FIG. 16: The power spectra of effective temperature fluctua- 
tion 6o + V&s and gravitational potential $ s — \]/ s at z = 1000 
in the Einstein and the Brans-Dicke model with SCDM pa- 
rameters. Normalization is arbitrary. 
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FIG. 17: The power spectra of effective temperature fluctua- 
tion Oo + and gravitational potential $ s — ^ s at z = 1000 
in the Einstein and the attractor model with SCDM param- 
eters. Normalization is arbitrary. 
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FIG. 18: The large scale power spectra of photon rms tem- 
perature fluctuation at the present in the Einstein and scalar- 
tensor models with SCDM parameters. Normalization is ar- 
bitrary. The scalar induced potential decay makes an uplift 
at the foot of the first Doppler peak and a bump around the 
present horizon scale. 



present horizon scales is higher than that in the Einstein 
model by uj^ 1 . However, the amplitude deviation inside 
the horizon (k ~ 1CT 3 Mpc -1 ) is somewhat reduced by 
this ISW (integrated Sachs- Wolfe) effect. On the other 
hand, the temperature fluctuation at the foot of the first 
Doppler peak which is blue shifted at the decoupling time 
is enhanced by the effect. This uplift is not so clear on 
the figure due to the smaller sound horizon length at the 
decoupling. In the attractor model, this effect is not so 
large on Sachs- Wolfe scales as on the scales around the 
first Doppler peak because the decay width is propor- 
tional to Ct>(</>) -1 at the horizon entry. 



V. OBSERVATIONAL QUANTITIES 

Observations of the large scale structure of the Uni- 
verse and the CMB anisotropies are the strong tools to 
constrain the theories of primordial fluctuation evolu- 
tion. Adopting scale invariant initial spectra, we present 
the matter power spectra and the CMB temperature 
anisotropy spectra based on the numerical calculations 
in the previous section. 

Before proceeding to the resulting spectra, we com- 
ment on the normalization. We find that, if the initial 
curvature perturbations are normalized to the same am- 
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plitude, the temperature anisotropy power on COBE nor- 
malization angular scales is almost model independent. 
However, this does not necessarily imply that the large 
scale power on the matter power spectrum is also model 
independent. The conservation of the adiabatic relation 
(A cr /A cm = 4/3) during when the mode is before the 
horizon entry is model independent. Here, the subscript 
r and to means a radiation and a non-relativistic com- 
ponent, respectively. According to Appendix, this fact 
indicates that, in the Brans-Dicke model, the familiar 
Sachs- Wolfe relation is modified as follows: 

6o + * s = i(l + |j)* s + 0(^ 2 ). (73) 

This is just an example. The actual modification also 
includes the contribution of the ISW effect and others. 



A. Matter power spectrum 

On subhorizon scales in the radiation-dominated 
epoch, the growth of density perturbations is inhibited 
by radiation pressure. After the transition from radia- 
tion domination to matter domination, perturbations on 
all length scales can grow by gravity. Hence in the mat- 
ter perturbation spectrum, there exists the imprint of the 
horizon scale at the matter-radiation equality time. We 
can see the difference in the horizon scales at the equal- 
ity time in Fig.[l] and Fig.^, The matter power spectra 
in the Einstein model and the scalar-tensor models are 
shown in Fig.[l9| and Fig.^o). It has been pointed out 
(in the context of Brans-Dicke cosmology pj|) that the 
variation of results in the larger expansion rate and 
then the turnover scale of the present-day matter power 
spectrum shifts to the smaller scale, accordingly there is 
more small scale power. It is found from the figures that 
large scale power deviates from its Einstein counterpart 
in proportion to Wq 1 - 



B. CMB temperature anisotropy spectrum 

The COBE normalized CMB anisotropy spectra are 
compared in Fig.|l] and Fig.^3|. Fig.^2| and Fig.|24| show 
their component decomposition which is calculated by 
directly projecting the inhomogeneity spectra at the de- 
coupling time into I space. 

First, we mention the shift of characteristic angular 
scales. The acoustic peak location (sound horizon scale) 
and the diffusion cut off scale are dependent on the hori- 
zon length at the decoupling time. Since the matching 
condition of <po constricts the deviation of the present 
horizon length, these angular scales directly represent the 
horizon length at the decoupling. Therefore, the shift of 
peak locations to smaller angular scales is due to the large 
gravitational constant at the decoupling. Although the 
damping scale also depends on the horizon length, it is 
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FIG. 19: The COBE normalized matter power spectra in the 
Brans-Dicke models with SCDM parameters. The turnover 
scale shifts according to the horizon scale at the matter radi- 
ation equality time. The difference in u) also affects the large 
scale amplitude. 
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FIG. 20: The COBE normalized matter power spectra in the 
attractor models with SCDM parameters. The turnover scale 
shifts according to the horizon scale at the matter radiation 
equality time. The large scale amplitude is almost the same. 



proportional to Ijj 5 . Hence the shift of damping tail lo- 
cation is not so large as that of peak locations. It means 
also that the width between the sound horizon scale and 
the diffusion cut off scale becomes thinner. 

Next, let us proceed to the variation of fluctuation am- 
plitude which arises from several effects. On the largest 
angular scale, the amplitude deviates by ujq 1 as in the 
matter power spectrum. The deviation on a Sachs- Wolfe 
scale is determined by the value of w(</>) -1 at its horizon 
entry In fact, the large scale tail in FigJ|^ deviates not 
so large as that in Fig.|24| which represents the value of 
o;(<^) -1 just at the decoupling time. Moreover, the scalar 
induced ISW effect is not negligible on large scales. We 
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FIG. 21: The COBE normalized CMB temperature 
anisotropy spectra in the Brans-Dicke models with SCDM 
parameters. 



FIG. 22: The component decomposition of the CMB tem- 
perature anisotropy spectra in the Einstein and the Brans- 
Dicke model which is calculated by directly projecting the 
inhomogeneity spectra at the decoupling time into I space. 
Normalization is arbitrary. Their curvature perturbations are 
initially normalized to be the same. 



can see in Fig.|25| that it damps the fluctuation amplitude 
on the scales which are red shifted at the decoupling time. 
The smaller sound horizon scale at the decoupling time 
also causes the attenuation of large scale amplitude in I 
space. As just previously noted, the amplitude on COBE 
normalization scales is almost model independent. Prob- 
ably the contributions of several effects cancel each other 
and the amplitude deviation around these scales is con- 
stricted. The foot before the first Doppler peak swells 
relatively to that in the Einstein model. The scalar in- 
duced decay contributes to this enhancement uplifting 
the originally blue shifted fluctuation, which can also be 
seen in Fig.^5| This decay boosts the acoustic peaks and 
the first peak is enhanced. Due to the peak location 
shift, each small scale peak enters the diffusion envelope 
more deeply than the corresponding peak in the Einstein 
model and thus it is damped more effectively. As an 
aside, comparing Fig. 2 and Fig. 23, we expect that the 
difference between the gravitational constant at the de- 
coupling G(z — 1000) and at the present time Go may 
be constrained as \G(z = 1000) - G \/G < 0.2. The de- 
tailed discussion of observational constraints of the model 
parameters is our future work. 



In the precedent study devoted to the CMB anisotropy 

the Brans-Dicke models 



in Brans-Dicke cosmology 20 
have higher acoustic peaks. On the contrary, according 
to our analysis, the Brans-Dicke models have generally 
lower acoustic peaks. In Fig]2| although the first peak 
height of the monopole component is almost the same as 
that in the Einstein model, the attenuation of the dipole 
component on that angular scale reduces the total peak 
amplitude. 
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FIG. 23: The COBE normalized CMB temperature 
anisotropy spectra in the attractor models with SCDM pa- 
rameters. 



VI. CONCLUSION 

We have comprehensively studied the perturbation 
evolution in the scalar-tensor cosmological model. It is 
shown that the scalar-tensor gravity requires some mod- 
ification to the standard evolution that is mainly caused 
by the variable gravitational constant and also brings ad- 
ditional process derived from its fluctuation. We provide 
the insight into their influence on matter evolution and 
thus it becomes possible to interpret the variation in the 
observable quantities. 

In the radiation-dominated epoch, the scalar-tensor 
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FIG. 24: The component decomposition of the CMB tem- 
perature anisotropy spectra in the Einstein and the attractor 
model which is calculated by directly projecting the inhomo- 
geneity spectra at the decoupling time into / space. Normal- 
ization is arbitrary. Their curvature perturbations are ini- 
tially normalized to be the same. 
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FIG. 25: The contribution of the scalar induced ISW effect to 
the CMB anisotropy spectra. The spectra not including the 
ISW contribution are calculated in the same way as in Fig.^j] 
and Figj24| and then the ISW contribution is added to them. 
The early ISW effect which is caused by radiation pressure is 
not included in any curve. 



gravity behaves like Einstein gravity. The influence on 
the perturbation spectrum is the larger horizon length 
alone. Its observational consequence is the smaller 
turnover scale in the matter power spectrum. On the 
other hand, in the matter-dominated epoch, the nature 
of gravity is qualitatively different from Einstein gravity. 
Consequently, aside from the shift of the acoustic peak 
locations and the diffusion cut off scale, the enhancement 
and the attenuation which originate from the scalar field 



perturbation and others appear in the CMB anisotropy 
spectrum. The attracting scalar-tensor gravity especially 
affects only the small scale power because of its attracting 
nature and the degree is far beyond the post-Newtonian 
deviation at the present epoch. 

In this paper, we have concentrated on studying the 
influence of the scalar-tensor gravity to cosmological per- 
turbations and their resultant observable power spectra. 
In the future work, we will put constraint on the model 
parameters ao and f3 by using the precise CMB and large 
scale structure data provided in near future. 
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APPENDIX: ANALYTIC APPROACH 

In this Appendix, the analytic approximate solutions 
of both the background and the perturbation equations 
in Brans-Dicke cosmology are presented (hence lo is con- 
stant here). Although these solutions are applicable only 
to Brans-Dicke models, we can use them for interpret- 
ing our numerical results of the variable lo models. The 
background cosmology is a flat SCDM model. To derive 
the approximate solutions, we expand the solutions per- 
turbatively by lo^ 1 . To the 0th order in lo^ 1 , we regard <f> 
as constant, = 1. (Therefore, X c is vanished to the 0th 
order in co~\) This means we fix the homogeneous solu- 
tion of Eq.(|12|) to the 0th order in lj~ x . Hence the 0th 
order equations and their solutions are the same as those 
in the Einstein model. The purpose of this section is to 
derive the corrections linear order in lo^ 1 . Higher order 
corrections are neglected here. (Hence this treatment is 
applicable only in large u> models) 

1. Background 

Several authors have studied cosmological solutions 
in Brans-Dicke cosmology [H- p5| . Here, we show the 
cosmological solution which is valid from the radiation- 
dominated epoch to the matter-dominated epoch. In the 
background cosmology, the linear order inw" 1 is assigned 
to (/>' , because the 0th order of (f> is constant and the force 
term of Eq.([l2|) is the linear order in to -1 . Since Eq.(|T2]) 
is the linear order in we can substitute 0th order 

solutions to the coefficient of each term and obtain the 
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linear order correction of the constant Oth order solution 
of (j>. In both the radiation-dominated epoch and the 
matter-dominated epoch, the homogeneous solution of 
Eq.(|l2"|) consists of <fi' = 0, which corresponds to another 
choice of G in Einstein gravity, and a decaying mode. 
More generally, we must give the boundary conditions for 
Hubble parameter and matter energy density to fix the 
magnitude of the decaying mode. However, we assume 
that the decaying mode is already negligible at z ~ 10 8 
(the initial time of our numerical calculations), otherwise 
it significantly affects the thermal history in the earlier 
epochs. (If the homogeneous decaying mode is dominant 
in (/)' , the energy density of <j> is proportional to a~ 6 . It 
is stiffer than that of radiation fluid.) So we employ, 
as the main contribution to the solution of Eq.([l2l), its 
particular solution that is 



, _ \/2 H eq a? eq 



1- I 1 + 

Q>eq 



l - 



2a 



,(A.l) 



eg , 



where a eq and H eq are the scale factor and Hubble param- 
eter without 0{ui~ l ) corrections at the matter-radiation 
equality time. Specifically, this choice determines the de- 
viation of Slo from unity. Integrating the above equation, 
we obtain 



4> = 0o + f{a) - /(a ), 



(A.2) 



where 



f(a) = - 



eq 



1 • •; 



(1 + 



2 \a eq + 



1 



(A.3) 



and an is the present scale factor. The boundary value 
O is 1 + _L + 0(w -2). 

a. Radiation- dominated epoch 



In the radiation-dominated epoch, Eq. (|A.3|) is reduced 



to 



/(«) 



In 4 




(A.4) 



We neglect 0(a 2 /a 2 q ) corrections here. In the period 
that a/a eq <C 1 is valid (i.e. early times in the radiation- 
dominated epoch), the last term can be neglected. This 
behavior of cf> is well known jll|. The contributions of 
(j)' and time variation of </> to a' /a are still smaller order 
than that of CDM energy density. This is because cf> is 
driven via the nonminimal coupling whose amplitude is 
regulated by the energy density of non relativistic matter. 
Then, we obtain the same relation between r\ and the 
expansion rate as that in Einstein gravity, 



V a 



(A.5) 



b. Matter-dominated epoch 



The solution in the matter-dominated epoch is, 
1 



f(a) = - \n(a/a eq ). 



So Eq.( |A.2[) becomes 



4> = 0o H ln(a/a ). 



The widely known power-law solution 

= 0o(^ 
Van/ 



(A.6) 



(A.7) 



(A. 



is reduc ed t o Eq.(A.7) with 0(ui 2 ) corrections. Using 
the Eq.(A/7), we obtain the following relation, 

2. Perturbation 

In this section, introducing a new parameter 

x = krj, (A.10) 

we derive the solutions of the perturbation equations. 
Obviously x is smaller than unity on superhorizon scales 
and larger on subhorizon scales. C' r and C m in the fol- 
lowing sections are proportionality constants. 

a. Radiation-dominated epoch 

We consider only the adiabatic initial conditions for 
matter perturbations in which T = 0. Here we also drop 
the anisotropic pressure correction to simplify the prob- 
lem. The Oth order initial conditions for matter pertur- 
bations are the same as those in Einstein gravity. Em- 
ploying the parameter x and neglecting the subdominant 
terms, Eq.(B9) is reduced to 



dx 2 



-— x c 

x dx 



a ( d . 



&lox a eq V dx 

4a, 

x 



(A.ll) 



in the radiation-dominated epoch. The Oth order solution 
of A c is 

A(°VC r = (^sin(a;/V3) -9cos(a;/V3)). (A.12) 

Then, as in the background case, we employ only the 
particular solution of Eq.(A.ll) which is 



27V3 



16ujx a e 



sin(a;/v / 3) 



^cos x — cos(x/v / 3) N 



(A.13) 
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(The particular solution is the only growing mode in X c . 
If the homogeneous mode exists in X c , it forms the grav- 
itational potential perturbation decreasing with time.) 
The behavior of this solution is as follows 



(x/V3) 



Superhorizon 
Subhorizon (^'^ 



Before the horizon entry, driven by the nonminimal cou- 
pling, X c grows proportional to x 2 . After the horizon 
entry, it stops growing and beg ins to oscillate. The first 
term in the parenthesis in Eq.( A.lj ) (cos a;) whose fre- 
quency is different from other terms represents the prop- 
agation of </> fluctuation and damps as the universe ex- 
pands. The oscillation of the remaining terms is driven 
by the matter acoustic oscillation which affects X c via 
metric perturbations. The factor l/\/3 in their argu- 
ment represents the radiation sound velocity. It is clear 
that the contribution of the additional terms in Eq.(|56|) 
and ( |60| ) to the gravitational potential is at most order 
(a/a eg )w _1 which is still smaller order than that of CDM 
perturbation. Hence, at the initial time early enough, the 
initial conditions for matter perturbations are not mod- 
ified. The relevant numerical solution is illustrated in 
Fig.|. 



b. Matter- dominated epoch 

In the matter-dominated epoch, the terms associated 
with c 2 or w are neglected, and then Eq.(p9|) is reduced 
to 



d 2 4 d 2 d 

TT^c H T-X-c + X c — T^c 

dx z x dx ujx dx 



JL Ac .(A.15) 



The modified evolution equation for A c becomes 



dx 2 



2 d 

x dx 



6 



2 d 

>jjx dx 



-A, 



4 

UJX 



-^-A C + ^X C + --^X C ~^X C . (A.16) 

dx* x dx rr 



6 

V 2 



Here we use Eq.( A.15 ). Obviously the homogeneous 
gro wing mode solution (the Oth order in oj^ 1 ) of 
Eq.( A.lq) is proportional to x 2 . Let C m be the propor- 
tion ality constant of it. Then, the particular solution of 
Eq.( |A.15| ) is 



X c l C m — — (l — 
uj V 

The asymptotic forms are 



3 sin x 3 cos x 



x 



(A.17) 



■^x 2 Superhorizon 
— . Subhorizon 



Before the horizon entry, it grows as in the radiation- 
dominated epoch. After the horizon entry, the propaga- 
tion mode damps away and the constant mode supported 
by the nonminimal coupling remains. The constancy of 
the surviving term is due to the constancy of the gravita- 
tional potential. This c onverg ence of X c is il lustrated in 
Fig.^. Substituting Eq.( |A.17 ) into Eq.( |A.16D , we obtain 
the modified solution of A c as follows 



A c / C m — X 



={»( 



cos x sin x 



+10 -x 2 



(A.19) 



The asymptotic forms of the correction terms are 



kz xi Superhorizon ^ ^ 



Subhorizon 



This correction is negligible until the horizon entry, and 
after the horizon entry the amplitude of the density per- 
turbation is slightly suppressed. Then, the solution of the 
gravitational potential with O(oj~ l ) accuracy becomes 



($ s - tf s )/C m = 12 




- 26 WA.21) 



At the horizon crossing, transient decay is induced by 
the suppression of X c growth. However, before and af- 
ter the crossing time, the linear order correction of the 
gravitational potential is flat as in the case of the Oth 
order. 



22 



Superhorizon 



— . Subhorizon 



(A.22) 



Finally, we mention the proportionality constants. Com- 
paring the same scale modes in the two models ever since 
the radiation-dominated epoch, we obtain the different 
C m s even if initially we set the C r s to the same value. 
C m itself contains the 0(w _1 ) correction which repre- 
sents the deviation formed around the matter-radiation 
equality time. On superhorizon scales, the correction co- 
efficient of the proportionality constant is 



c bd i c e ^ 1 , 

° m /Cm ~ 1 + lOw 



(A.23) 



(A.18) according to our numerical calculations. 
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